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Abstract 


Various  problems  concerning  the  effects  of  the  boundaries 
of  the  ocean  on  the  propagation  of  pressure  waves  in  the  ocean 
are  considered.  The  propagation  of  a  transient  pressure  wave  in 
a  wedge  shaped  region  of  fluid  is  treated.  This  is  the  model 
chosen  to  describe  the  situation  in  which  an  underwater  explosion 
takes  place  in  a  coastal  ocean  region  which  is  characterized  by  a 
strongly  sloping  bottom.  In  an  attempt  to  study  the  effects  of  the 
polar  ice  cap  on  the  propagation  of  a  pressure  wave,  the  reflection 
of  a  plane  wave  onto  a  rough  boundary  separating  a  fluid  half  space 
and  a  thick  fluid  layer  of  differing  sound  speed  and  density  is  con¬ 
sidered.  These  results  are  currently  being  used  to  construct  the 
response  to  a  transient  pressure  pulse  and  to  generate  numerical  re¬ 
sults  for  conditions  representative  of  underwater  explosions.  The 
final  section  presents  numerical  values  of  the  reflection  coefficient 
as  a  function  of  grazing  angle  for  the  case  of  a  plane  wave  incident 
on  a  porous  elastic  bottom.  The  analytical  expressions  used  were 
derived  in  an  earlier  report  of  this  series. 
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I.  INTRODUCTION 


The  underwater  detonation  of  an  explosive  charge,  whether  it  be  nuclear  or 
conventional,  produces  a  region  of  high  pressure.  This  pressure  ia  transmitted  to 
the  surrounding  water  and  propagates  as  a  shockwave  in  all  directions.  It  is  well 
known  that  the  geometric  and  mechanical  properties  of  the  boundary  media  as  well 
as  the  propagating  media  itself  can  strongly  modify  the  pressure  pulse  that  is 
sensed  at  a  given  location.  The  study  of  such  effects  is  necessary  for  an  under¬ 
standing  of  underwater  explosion  effects  and  explosive  echo-ranging. 

The  actual  situation  is  of  course  quite  complex.  In  the  immediate  vicinity  of 
the  explosion  non-linear  effects  are  important  due  to  the  extremely  high  pressures 
generated.  The  pressure  pulse  propagates  through  water  which  is  Inhomogeneous, 
the  sound  speed  and  density  varying  with  depth.  The  boundary  media  can  themselves 
be  quite  conpllcated  in  their  mechanical  behavior  and  furthermore  the  interface 
between  the  water  and  the  boundary  is  often  highly  irregular.  In  order  to  make 
any  mathematical  analysis  tractable  many  simplifying  assumptions  must  be  made. 

For  the  purposes  of  this  report  we  restrict  ourselves  to  analyzing  the  pressure 
field  in  terms  of  linear  acoustic  theory,  i.e.,  at  distances  large  enough  from 
the  source  of  the  explosion  so  that  the  pressure  amplitude  has  decreased  suffi¬ 
ciently  to  permit  the  use  of  the  linearized  equations  of  motion.  The  fluid  through 
which  the  pressure  pulse  propagates  is  assumed  homogeneous. 

An  earlier  report^  concerned  Itself  mainly  with  the  effects  of  the  mechanical 
properties  of  the  boundary  mdela,  treating  the  boundary  media  as  a  liquid  satu¬ 
rated  porous  elastic  solid  (mathematical  model  of  a  consolidated  sedimentary 
bottom).  In  this  study  the  interface  between  the  bottom  media  and  the  water  was 
taken  as  a  smooth  plane  parallel  to  the  free  surface  of  the  ocean. 

The  present  report,  for  the  most  part,  will  emphasize  the  geometric  effects 
of  the  boundaries  rather  than  the  mechanical  behavior  of  the  boundary  media. 

Sect.  II  treats  the  case  of  a  pressure  pulse  propagating  in  a  coastal  region  char¬ 
acterized  by  a  strongly  sloping  bottom.  Sect.  Ill  concerns  itself  with  the  re¬ 
flection  of  plane  pressure  waves  from  a  rough  boundary  while  Sec.  IV  treats  the 
reflection  of  waves  from  sinusoidal  boundaries  between  various  types  of  media. 

These  results  are  currently  being  used  to  find  the  response  to  a  transient  pressure 
wave  and  to  generate  numerical  results  for  conditions  representative  of  underwater 
explosions.  Finally,  Sec.  V.  presents  numerical  results  for  the  reflection  of  a 
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plane  wave  incident  from  a  liquid  onto  a  liquid  saturated  porous  elastic  solid. 
(The  theory  of  this  was  presented  by  Blchler  and  Rattayya.1) 

II.  PRESSURE  FUIflE  PROPAGATION  FROM  AN  UNDERWATER  EXPLOSION  IN  THE 

VICINITY  OP  A  SLOPING  BOTTOM 

In  recent  years  interest  has  developed  in  the  propagation  of  pressure  waves 

from  underwater  explosions  taking  place  in  coastal  ocean  regions.  These  areas  are 

characterized  by  a  strongly  sloping  bottom  consisting  of  volcanic  or  coral  rock. 

The  angle  that  the  bottom  makes  with  the  horizontal  varies  from  about  10°  to  as 
o  2 

much  as  40  .  The  situation  that  these  regions  present  is  quite  different  from 
that  of  either  the  continental  shelves  or  the  deep  ocean  basins,  where  the  water 
depth  can  be  considered  to  be  fairly  constant. 

In  order  to  study  the  problem  of  underwater  explosions  in  the  vicinity  of 
strongly  sloping  bottoms,  the  situation  is  idealized  by  considering  a  transient 
point  or  line  source  situated  in  a  wedge  shaped  region  occupied  by  a  slightly 
conqpressible ,  non-viscous  liquid  of  density  p,  and  constant  sound  speed  c.  Cylin¬ 
drical  coordinates  (r,0,z)  are  introduced  such  that  the  z-axis  coincides  with  the 
edge  of  the  wedge  and  the  boundary  walls  of  the  wedge  are  given  by  6-0,  and  0-f. 
The  boundary  at  0-0  is  taken  to  represent  the  ocean  surface,  while  the  boundary 
at  0-f  will  correspond  to  the  ocean  bottom.  The  point  source  is  located  by  the 
cylindrical  coordinates  (r  ,0Q,o)  as  shown  in  Fig.  1.  In  contrast  to  the  problem 
of  a  point  source  located  in  an  ocean  of  constant  depth,  which  benefits  from  cy- 
lindricalsymmetry,  the  problem  of  a  point  source  in  a  wedge  shaped  region  is  truly 
a  three-dimensional  problem. 

A  rigorous  mathematical  model  for  the  above  situation  should  assume  that  the 
bottom  boundary  0*f  is  in  contact  with  an  elastic  solid  or  another  liquid  of 
different  density  and  sound  speed,  but  this  problem  has  been  shown  to  be  mathe¬ 
matically  Intractable.  K^arsley's  thesis  gives  a  detailed  description  of  the 
difficulties  involved.  Therefore  as  a  first  approximation  to  the  problem  the 
ocean  bottom  will  be  assumed  to  be  infinitely  rigid,  i.e.,  the  normal  derivative 
of  the  pressure  vanishes  at  0* f,  and  the  ocean  surface  will  be  taken  as  a  free 
surface,  i.e.,  the  pressure  is  zero  at  the  surface  0*0.  The  method  used  will  very 
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closely  follow  that  of  Oberhettinger.  This  reference  also  contains  an  excellent 
bibliography  of  earlier  work  on  the  problem  of  wave  diffraction  by  wedges.  The 
present  work  goes  beyond  Oberhettinger,  in  that  it  solves  the  problem  with  one 
boundary  rigid  and  one  boundary  free,  whereas  Oberhettinger  considers  the  case 
where  both  boundaries  are  either  free  or  rigid. 


A.  Line  8ource  Bx citation 


First  consider  the  case  of  a  time  dependent  line  source  parallel  to  the  z- 
dlrection  situated  on  the  line  S (rQ,0  )  within  the  wedge.  The  pressure  field 
within  the  wedge  0  <  0  <  f  must  then  satisfy  the  wave  equation 


is.  1  dg  .  1_  i%  i_&  -6(r-ro)t(a.9o)f(t) 

dr2  75?  r2fc2'c2at2‘ 


(n.i) 


where  f(t)  represents  the  time  history  of  the  source  and  6(x)  is  a  Dirac  delta 
function,  together  with  the  boundary  conditions 


p(r,e)  -  0 

l(r,0  -  0 


(II. 2) 


Taking  a  Laplace  transform  in  time  of  both  sides  of  Eq.  II. 1  results  in 


where 


and 


dr2  ra*2  c2’ 


t(r-ro)s(e^o)F(s) 


(H.3) 


OB 

P (r,e)  •  J  p(r,0,t)"stdt  \ 


Y*i« 


(II. 4) 


p(r,0,t)  -  gJI  J  p(r,0)est  ds 

4-i« 


p(s )  -  J  f(t)est  dt 


(II -5) 


is  the  Laplace  transform  of  the  time  dependent  source  strength.  In  order  to  solve 
Eq.  II. 3  together  witfr  the  transformed  boundary  conditions  corresponding  to 
Eq.  II. 2,  consider  first  the  following  problem;  a  transient  source  is  located  at 
r*rQ,0"$o  in  a  wedge  shaped  region  0  <  6  <  a,  the  boundaries  of  which  are  pressure- 
release  surfaces,  i.e.,  the  pressure  vanishes  at  0-0  and  at  0>a.  The  solution  to 
this  problem  can  be  easily  derived  (see  Appendix  A)  and  is  given  by 


p(r 


..)  -  .. 


la  /  •Il(«>»»n(e-8o)*co,”,n(9Wo)]I»  (t  r  <  >  K„  ^  >  ) 

W  n  n 


(II.6) 
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where  r  <  and  r  >  are,  respectively,  the  smaller  or  larger  of  the  quantities  r  and 

rQ,  and  «o=  1,  «n=  2  for  n  >  1.  The  solution  for  an  image  source  located 

at  r=r  6=cc-d  is 
o  o 


p(r  ,9 )  = 


2U1  ^ 

2a  / 
n=0 


en[cosun^'a+eo^co8un^‘0o4a)]1u  r  <)Kv  (c  r 

n  n 


(II.7) 


Combining  the  two  solutions  and  letting  a  ■  2$  one  arrives  at  the  Laplace  trans¬ 
form  of  the  solution  to  the  problem  posed  by  Eqs.  II. 1  and  II. 2,  and  this  can  be 
written  as 


00 


p(r,e)  =  €n^cosun^0“eo^  "  cosun(e+0o) 


(II. 8) 


+  cosun(Weo-2*)-coBi>n(e-eo+2t))Ii)  (|  r  <)KU  (f  r  >) 

n  n 


For  convenience  in  further  manipulation  define  the  function  G  by 


5(r,0;ro,0o) 


1 


Y.  *n[cOBUn(e'9o) 


n=0 


cosv  (0+0  )]I  (-  r  <)K  (-  r  >) 
nv  o  v  c  v  c 

n  n 


(II. 9) 


so  that 


p(r,0)  =  {G(r,0;ro,0o)  +  5(r ,0 ; ro,2rjf-0Q ) }  F(s) 
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Following  Oberhettinger  it  can  be  shown  that  G  can  be  put  in  the  form 
G(r,0;ro,0o)  =  L(r,ro,0-0Q)  -  £(r,ro,0+0Q) 

where 


(II. 10) 


(II. 11) 


h 


£(*,VP) 


c. 

kY  Kolt 


r  2^  2 
[r  +r 

o 


2  rrQ  cos(4rujr  +  tp)]*) 


n«n. 


2,  2 
r  +r 

o 


+  2 


rrQ  coshx)^]  X 


(11.12) 


I' 


sin  (n-cp) 


sln[|^(n+9)] 


^cosh(~)-cos[|^(jt-«p)]  cosh(^)-cos[|^(jt+9)) 


dx 


where  is  the  largest  positive  or  negative  number  which  is  smaller  than  or  equal 
to  -  (jt-Kp/4^)  and  is  the  largest  positive  or  negative  number  which  is  smaller 
than  or  equal  to  («-<p/4t).  If  is  greater  than  the  sum  in  Eq.  11.12  is  zero. 


B.  Point  Source  Excitation 


The  line  source  solution  derived  in  the  previous  section  can  be  used  to  derive 
the  solution  for  a  point  source.  Using  Cylindrical  coordinates  (r,0,z)  the  point 
source  is  located  at  (r^e^o ).  Eq.  II. 1  now  becomes 


dr2  rS?  r2*2  dz2  =2&t2 


6(r-ro)6(e-0o)6(z)f(t) 


(11.13) 


Taking  a  Fourier  Cosine  transform  in  z  and  a  Laplace  transform  in  time  gives 


+  1  dg  +  i  *2S  -«2  =>-  «(r-rJ4(e-ejF(.) 


dr* 


1  ^ 
r  dr 


r  dfi  c 


2r 


(11.14) 


where  \  is  the  Fourier  Cosine  transform  parameter.  It  is  then  obvious  that  one 

can  derive  the  point  source  solution  from  the  line  source  solution  by  replacing 

2  g  2  i 

s/c  by  (s  /c  +\  )2  in  the  line  source  solution,  multiplying  each  term  of  the  solu¬ 
tion  by  (l/«  cosXz)  and  integrating  with  respect  to  \  from  0  to  ».  In  performing 
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the  above  operation  one  can  use  the  following  well-known  relationship 


(r2+z2)~^exp[~  (r2+z2)2] 


+  \2)^]  cos\z  dX 


(11.15) 


5 


Operating  on  Eq.  11.12  in  the  above  described  manner  we  obtain 


<* 


1(0) 


■hf 


exp{-  ^[r2+r2+z2-2rrQcos(0+4nf )  ]*} 
[r  +r^+z  -2rrocos  (0+4ni|; )  ]  * 


l. 


n=n. 


'  H<x) 


sin[|^(n-9)3 


sin[|^(g->-9)3 


cosh(||)-cos  [|^(x-0  )  3  cosL(~)-cos  [§^(*+0  )  3 


-  dx 


where 


H(x) 


8  2  2  2  i 

exp[ - (r  +r  +z  +2rr  coshx)*3 

r  c  v  o  o 

pop  i 

(r  +r  +z  +2rr  coshx)^ 

'  o  o 


Finally  using  Eqs.  II. 10  and  II.  11  p(r,0)  can  be  written  in  the  form 


(11.16) 


(11.17) 


p(r,0)  -  [t(r,ro,0-0o)-L(r,ro,e+0o)  +  L(r ,ro,0-2t+0o)-L(r,ro,0+21Jr-0o) 3  F(s)  (II. 18) 

where  L(r,rQ,0)  is  defined  by  Eq.  II.  16.  Each  L(r,rQ,0)  consists  of  two  kinds 
of  terms,  one  a  sum  of  terms  each  of  which  has  the  same  functional  dependence  as 
the  point  source  but  located  at  image  points  of  the  source  with  respect  to  the 
boundaries  of  the  wedge.  The  number  of  images  depends  on  the  numbers  n^  and  n^ . 

The  other  term  is  given  by  an  integral  expression  which  vanishes  whenever  the 
angle  of  the  wedge  $  is  a  submultiple  of  n,  i.e.,  ♦  ■  x/m,  m  =  1,2....  If  this 
is  the  case  the  solution  is  then  given  by  the  effect  of  the  original  source  and  a 
finite  number  of  image  sources  which  are  necessary  to  satisfy  the  specific  boundary 
conditions.  The  integral  expression  has  been  called  the  edge  diffracted  wave  by 
Biot  and  Tolstoy. 


C.  Transient  Solution 

It  has  been  assumed  in  writing  Eqs.  II. 1  and  11.13  that  both  the  line  and  point 
source  are  time -dependent,  and  that  their  time-dependence  is  given  by  f(t).  In  ao- 
plying  the  Laplace  transform  to  the  function  f(t)  it  has  been  assumed  that  f(t) 
vanishes  for  t  <  0.  What  remains  to  be  shown  is  how  the  solutions  depend  on  f(t). 

In  order  to  present  results  for  the  time -dependent  source  it  is  convenient  to 
introduce  the  following  notation 

Rn(0)  =  {r2+r2+z2-2  rrQ  cos(0+4nt))^  (II.19) 
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1 


ppp  ^ 

R  (x)  =  {r  +r  +z  +2  rr  cosh  x}4 
nx  '  *•  o  o 


(11.20) 


g(x,0) 


sin[|^(«-0)j 


sin[f-(x+e)] 
+  - 2* - 


cosh  (||)  -  cos  [|^(«  -0 )  ]  cosh(j^)  -  cos  [|r(x+0  )  ] 


Using  the  above  notation  I(r,rQ,0)  can  be  written  as 


(11.21) 


£ 

i  V”  -  s/c  R  (0)  .  ‘ 

t(r,r0,S)  -  TJJ  R  CeJ - I 

n«n. 


-  s/c  R  fx) 


6  Rn(x) 


g(x,0)  dx 


(11.22) 


The  solution  corresponding  to  the  case  F(s)«l,  where  F(s)  as  given  by  Eq.  II. 5 
is  the  Laplace  transform  of  the  function  f(t),  represents  the  response  of  the 
system  to  a  Dirac  delta  excitation  applied  at  the  time  t=0. 

This  response  is  given  by 


"2  R 

,  V  6[t-  -S 

P4(r>s.t) 

n«l 


6[t. 


c 

O' 


r  (e-2t+e  )  r  (e+2*-e  ) 

e[t-  -a— — — j  6[t.  _s — ^ — — ] 


Rn(0-2*+0o) 


Bn(e+at-V 


(n.23) 


where 


R  (x) 

X1  6[t — — ] 

'  l^tn  l  - - ^(x,0-0o)-g(x,0+0o)+g(x,0-2t+0o)-g(x,0+2tlr-0o))  dx 


n' 


2.2  2  2  2 
c  t  -r  -r  -z 

xx  »  cosh"  ( - 2rr — - )  • 


(11.24) 


(11.25) 


For  any  given  f(t)  the  first  four  summations  can  be  obtained  explicitly  while 
the  integral  contribution  can  be  evaluated  numerically. 


7 


It 


III.  SCATTERING  OF  A  PLANE  WAVE  BY  A  ROUGH  NON -SINUSOIDAL  BOUNDARY  BETWEEN  A 

SEMI -INFINITE  FLUID  MEDIUM  AND  A  THICK  FLUID  LAYER 

The  problem  considered  in  this  section  is  that  of  a  plane  wave  scattered  by 
a  rough  non-sinusoidal  boundary  between  a  semi-infinite  fluid  and  a  thick  fluid 
layer.  This  is  the  first  step  necessary  for  the  consideration  of  the  scattering 
of  a  transient  pressure  pulse  by  the  Arctic  icecap.  These  results  are  currently 
being  used  to  construct  the  response  to  a  transient  pressure  pulse  and  to  generate 
results  for  conditions  representative  of  underwater  explosions . 

The  reflection  of  waves  from  a  periodic  boundary  has  been  studied  quite  ex¬ 
tensively  and  for  a  discussion  of  this  problem  see  the  following  section.  Reflec¬ 
tion  of  waves  from  arbitrary  non-sinusoidal  boundaries  has  not  been  studied  so 
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extensively.  Recently  Abubakar  and  Dunkin  and  Eringen  have  studied  the  reflec¬ 
tion  of  waves  from  an  arbitrary  non-sinusoidal  surface  by  means  of  a  perturbation 
method.  It  is  essentially  this  method  which  is  used  in  the  present  analysis. 

Consider  the  case  of  a  plane  wave,  propagating  through  a  fluid  of  density  pQ  and 
sound  speed  cq,  incident  on  a  liquid  layer  of  density  p^,  sound  speed  c^  and  of 
mean  height  h.  The  interface  between  the  fluid  and  the  liquid  layer  is  taxen  to 
be  z-*f(x)  where  c  is  a  small  parameter  and  f(x)  is  the  arbitrary  function  of  x. 

The  fluid  lies  in  the  half  space  z  >  cf(x),  while  the  liquid  layer,  whose  upper 
boundary  (z»  -h)  is  taken  to  be  planar,  then  lies  in  -h  <  z  <  ef(x)  (Fig.  2),  where 
z  is  measured  positive  into  the  fluid  half  space. 

As  before,  a  plane  sound  wave  of  unit  amplitude  is  «RRnmf»d  to  be  obliquely 
incident  upon  the  irregular  boundary  separating  the  two  fluids.  The  velocity 
potential  of  this  wave  is  given  by 

<t±  -  e1(Co*-V>  (III.l) 

where  ao=kQsin0o,  po=kQ  cos (kQ=(o/co),  the  projections 
number  on  the  x-  and  z-axis  respectively,  is  the  angle 
mal  makes  with  the  z-axis  as  shown  in  Fig.  2. 

The  reflected  wave  <Pr  is  assumed  to  be  of  the  form 

<Pr  -  «  J  A(?) 

-00 

v  * — 2f  2< 

where  yQ  =  '  k0"5  •  *n  the  case  of  a  sinusoidal  boundary 
cites  a  denumerably  infinite  set  of  plane  reflected  waves 
is  assumed  to  scatter  a  continuous  spectrum  of  plane  waves. 


of  the  incident  wave 
that  the  plane  wave  nor- 


d?  (III. 2) 

the  incident  wave  ex- 
whereas  in  this  case  it 
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The  transmitted  wave  Cp.  is  assumed  to  be  of  a  similar  form 

X0 


and  is  given  by 

(HI. 3) 


+  i_  f  B(!)[e1(!x'yl2)  +  D  el(Jx+Ylz)]  d? 

*QO 


where  0^  is  the  angle  of  refraction  of  the  transmitted  wave,  a^«k^  sinS^, 
H-^k-cosA.  are  the  projections  of  the  transmitted  wave  number  on  the  x,  and  z  axes 

1  -.7  2  2” 

respectively,  and 

We  assume  that  the  upper  surface  of  the  liquid  layer  z«  -h  is  a  pressure- 
release  surface,  i.e.,  the  potential  given  by  Eq.  III. 3  must  vanish  at  z=  -h. 
Using  this  condition  we  can  solve  for  C  and  D  in  Eq.  III. 3  and  we  obtain 


C  -  -  e2i“lh 
D  -  -  e21Vl» 

Substituting  these  values  back  into  (III. 3)  ve  obtain 
<pt  -  T  eUo1x-u1z)[l  e21li1(Z+h)]  + 

B(S)  ei(«*-Vi*)[1.e2iV1(*+h)]  dJ 


(III.4) 


(III.?) 


Hie  rough  boundary  separating  the  two  fluids  is  given  by 
z  »  cf(x) 


(III. 6) 


where  e  is  a  small  characteristic  parameter  (for  example  the  ratio  of  the  maxi¬ 
mum  amplitude  of  the  boundary  to  the  smaller  of  the  two  acoustic  wavelengths 


2itc 


2nc, 


’(D  1  0) 


or  X,  «  -rr*). 

In  order  to  determine  the  unknowns  R,T,A(?),  B(?)  appearing  in  Eq.  III. 2  and 


Eq.  III. 5  the  boundary  conditions  at  the  interface  z«*f(x)  must  be  satisfied.  The 
two  boundary  conditions  are  continuity  of  pressure  and  continuity  of  normal 
velocity.  In  terms  of  the  potential  functions  the  condition  of  continuity  of 
pressure  is 

p0(<Pi  +  <Pr)  -  px9t  J  Z  -  ,f(x)  (in. 7) 
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•  2  -4 

The  components  of  the  unit  normal  to  the  surface  z=*f(x)  are  n  =  -cf'[l+(ef') 

2  — —  ^ 

n  “[l+Ccf7)  ]  2,  where  f'=df/dx,  and  so  the  normal  velocity  condition  is  written  es 


d(c?.-Kp  )  dcp  dtp  dip 

i  r  _/  /  l  .  r\  t 


-  *f' (sr  +  sr)  -  sr  -  *f '  sr  *  *  =  «f<x> 

As  mentioned  earlier  «  is  taken  to  be  a  small  parameter  and  therefore  we 
expand  the  functions  A($)  and  B(j)  into  power  series  in  t* 


(III. 8) 


OD 

A(5)  c\(S)  -  (5)  +  ,2Ag(5)  +  ... 

n«l 


(III. 9) 


B(5)  -  V  c\(?)  -  cB^!)  +  €%({)  +  ... 
n*l 


(III. 10) 


Substituting  these  expressions  into  Eq.  III. 2  and  Eq.  III. 5  and  then  the  resulting 
expressions  into  Eq.  III. 7  yields 

« 

Y.  ‘Vl(v°,f+Sx)  d:}  ■ 


p  (Te1(ffl3t-,1ltf)[l-e2i,1l(€f+h)]+  f-  f  V  ‘"B  el(°lX‘VlEf)[l^2lYl(,f+h)J  dg } 

11  -eo  n 

n«l 


(III. 11) 


Now  expanding  the  exponentials  of  the  form  e0***  into  a  power  series 


“«f . !  +  ssi  +  i2£.  + 

1  *  1!  +  2!  +  "• 


(III. 12) 


and  equating  the  coefficients  of  like  powers  of  c  in  Eq.  III. 11  we  obtain  the 
equations 

ia  x 


PoReiaox  -  PlT(i.e21^lh)  el0lx  =  .  pQe10o 


(III. 13a) 


^  J  A1(?)el5x  d;  -  J  B(g)  (l-e2lYlh)  e15x 

—ao  —oo 


<*5 


i  f[pouo(l-R)  eiao*-  -  p1M1T(l+e2i^lh)  eialX]  ,  etc. 


(III. 13b) 


Eqs.  III. 13a  and  III. 13b  are  determined  from  the  zeroth  and  first  order  e  terms 
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J 


In  Eq.  III.  11,  respectively.  In  order  that  Eq.  III.  13a  be  satisfied  for  all  values 
of  x  it  is  necessary  that  which  means  that 


sind 

o 


c 

o 


8in01 

C1 


(III. 14) 


which  is  the  familiar  expression  of  Snell's  Law. 

In  the  same  manner  one  satisfies  the  normal  velocity  condition  Eq.  III. 8  and 
obtains  the  following  two  equations  corresponding  to  the  zeroth  and  first  order 
terms  in  <• 

R  +  1^(1  +  e*a^x  T  «  noe*aox  (III. 15a) 


If  I  Vl*15*  +  If  /  V1U+«2i^lh)  B  eltx  d5 

■08  *00 

(111.15b) 

-  (aQf/-i^f)(l+R)  eiCTox  -  (a1f/-i^f)(l-e2iulh)  TeialX 
Making  use  of  Eq.  III. 14,  Eqs.  III. 13a  and  III. 13b  become 

PoR  -  p1(l-e21,1lh)  T  -  p0  (III. 16) 

PoR  +  M1(l+e21,1lh)  T  -  po  (III.17) 

Solving  for  R  and  T  we  obtain 


£  (l-e21tllh)  -  £  (l+e21wlh) 

R  -  — - — -  (III. 18) 

—  (l-e2i^lh)  +  —  (i+e2i^lh) 

^1  4o 


-i 

T  -  — - - -  (III. 19) 

—  (l-e2i^lh)  +  --  (l+e2i4lh) 

*1  ^0 

In  order  to  solve  for  the  Amotions  A(?)  and  B(§),  multiply  ~ach  of  the 

•  if'  jr 

Eqs.  III.13b  and  III. 15b  by  e  and  then  integrate  over  x  from  -®  to  ®,  making 
use  of  the  identity 

J  e1^"5  ^x  dx  -  2*6(5-$')  (III. 20) 
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where  6(0  is  the  Dirac  delta  function.  The  transformed  Eqs.  III. 13b  and  III. 15b 
become,  respectively 


Pi(l-e2lYlh)  B1(S)  - 

(III. 21) 

1  f(5-o0)(p0Mo(l-R)  -  p1n1T(l+e21^h)} 


where 


voA1(?>+  V1(l+e2lYlh)  - 

(III. 22) 

lf(5-Po)t(ao(!-00)-P2](l+H)-[<Jo(5-oo)-p2Kl-e2lMlh)  T) 


Hi)  -  I  f(*)  e‘i{X  <bc 

_eo 

(hi. 23) 

f(*)  •  gj  J  Hi)  e1!x  IS 


Eqs.  III. 21  and  IU.22  can  be  solved  for  and  B^.  The  present  analysis  will 
only  calculate  the  first  order  perturbation  to  the  reflected  wave  and  so  we  need 
only  obtain  which  is  given  by 


Ajft)  -  21  r(t-°  HU-  jv2  tan  n^h][l- 

poMl 


ip  V  -i 

-2-2  tan  y^]}"1  X 

poYl 


(III. 24) 


Pt  Pt  t  Pi  o  Pi  o  * 

(n  (l - -)  + - tan  y.h  tan  n,h[o  (5 -a  )(l - )-pf+  ]) 

ov  p  '  p  p.y,  T1  T.  ov*  o,v  p  7  p  toJJ 

Ko  o  1T1  Ko  Ko 


Insertion  of  this  value  of  A1  into  Eq.  III. 2  yields  the  first  order  pertur¬ 
bation  to  the  reflected  potential  as 

<Prl)  "  h  J  *1(5)  ei(V+!x)  a?  (III.25) 

•00 

V~2  21 

k  -5  .  This  integral  can  be  evaluated  asymptotically,  i.e.,  at  dis- 

00  2  2 
tances  from  the  origin  that  are  large  enough  so  that  kQr  »  1,  where  r*ic  +z  and 

the  result  is  (see  Appendix  A) 

^r^  ~2*r  el^°r"  ^  sin0)  0080  *  (III. 26) 
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From  Eq.  III. 2  it  is  clear  that  the  first  order  perturbation  to  the  re¬ 
flected  potential  is  given  by  a  cylindrically  spreading  wave. 

IV.  REFACTION  OF  PLANE  WAVES  AT  AN  IKREGUIAR  BOUNDARY 

The  reflection  of  a  plane  sound  wave  from  a  non-planar  interface  between  two 

8 

media  has  been  termed  by  Uretsky  a  "marvelously  complex  problem. "  The  problem 
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was  first  attacked  by  Rayleigh  and  has  since  apparently  captured  the  fancy  of 
many  workers  for  there  is  a  wealth  of  literature  dealing  with  the  problem.  The 
bulk  of  this  literature  however  is  concerned  with  the  simplest  case  which  is  that 
of  a  sinusoidal  pressure -release  boundary.  This  simplest  case  is  trouble  enough 
however,  for  its  solution  involves  solving  infinite  sets  of  equations  for  which, 
of  course,  only  approximate  solutions  have  been  obtained.  There  is  much  less  lit¬ 
erature  dealing  with  the  case  of  the  non-perfectly  reflecting  boundary  such  as  a 
sinusoidal  boundary  between  two  fluids  of  different  characteristic  impedance,  no 
doubt  because  the  problem  is  a  degree  more  difficult  than  the  perfectly  reflecting 
boundary  inasmuch  as  the  amplitudes  of  the  rays  transmitted  into  the  second  fluid 
constitute  a  second  infinite  set  of  unknowns  to  be  dealt  with. 

We  will  briefly  mention  the  other  analyses  of  reflection  from  a  non-perfectly 
reflecting  sinusoidal  boundary.  Miles10  treats  the  case  of  oblique  incidence  upon 
a  sinusoidal  boundary  between  two  semi-  infinite  fluid  media  in  contact  and  both  he 
and  Grasyuk11  also  look  at  the  case  of  oblique  incidence  upon  the  sinusoidal  bound¬ 
ary  between  a  semi-infinite  fluid  in  contact  with  a  semi-infinite  elastic  medium. 

In  both  these  analyses,  the  amplitude  of  the  undulations  of  the  boundary  is  assumed 

to  be  small  compared  to  a  wavelength  of  sound  while  the  slope  of  the  boundary  is 
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considered  to  be  non-negligible .  However  it  has  been  pointed  out  by  Lapin  that 
Grasyuk' s  work  has  errors.  Rayleigh  in  his  work  took  the  opposite  point  of  view, 
i.e.,  the  slope  of  the  boundary  is  assumed  small  while  the  amplitude  of  the  undu¬ 
lations  is  not.  Rayleigh's  analysis  is  also  available  in  Ref.  13.  Rayleigh  lim¬ 
ited  his  attention  to  normal  incidence. 

In  the  analysis  that  follows,  we  extend  Rayleigh's  work  to  allow  for  oblique 
incidence.  An  approximate  solution,  entirely  analogous  to  that  of  Rayleigh,  is 
obtained  under  the  conditions  that  the  acoustic  wavelength  is  much  less  than  the 
wavelength  of  the  undulations  of  the  boundary  and  the  slope  of  the  boundary  is 
negligible.  We  then  treat  in  turn  the  problems  of  reflection  from  a  sinusoidal 
boundary  of  a  plane,  obliquely  incident  wave  between  l)  two  semi-infinite  fluid 
media,  2)  a  semi-infinite  fluid  and  a  semi-infinite  elastic  medium,  3)  a  finite 
thickness  fluid  layer  sandwiched  between  two  different  semi-infinite  fluid 
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media,  and  4)  a  finite  thickness  elastic  layer  sandwiched  between  two  semi- 
Infinite  fluid  media.  For  these  four  cases,  an  approximate  solution  Is  again  ob¬ 
tained  under  the  restriction  Just  mentioned,  and,  for  the  last  two  cases,  with  the 
additional  restriction  that  the  thickness  of  the  layer  be  much  greater  than  the 
amplitude  of  the  boundary  undulations. 

Before  beginning  the  analysis,  a  word  must  be  said  about  the  controversy  that 

exists  regarding  Rayleigh's  assumption  that  the  reflected  and  transmitted  waves 

can  be  written  as  a  superposition  of  plane  waves,  each  propagating  In  a  discrete 

direction.  Although  it  is  agreed  that  this  Is  a  valid  representation  for  the 

scattered  field  at  a  large  distance  from  the  boundary,  the  controversy  pertains 

to  the  question  of  whether  this  representation  is  valid  near  the  boundary,  that  is, 

in  the  region  between  the  extremes  of  the  undulations.  No  rigorous  proof  that 
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Rayleigh' 8  assumption  is  incorrect  has  yet  been  published  although  Uretuky  pre- 

15 

sents  a  strong  argument  to  that  effect.  Marsh,  on  the  other  hand,  has  published 
an  argument  which  he  claims  Justifies  the  Rayleigh  assumption.  In  the  analysis 
that  follows ,  we  unashamedly  and  without  providing  any  additional  Justification 
use  the  Rayleigh  approach. 

A.  Scattering  from  A  Sinusoidal  Boundary  between  Two  Semi-Infinite  Fluid  Media 

The  geometry  of  the  problem  is  shown  in  Fig.  3*  The  notation  used  closely 
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follows  that  of  Rayleigh  or  Orasyuk.  A  plane  sound  wave  of  unit  amplitude* 
(defined  by  the  velocity  potential  <p^)  is  assumed  to  be  obliquely  incident  upon 
the  sinusoidally  rough  boundary  (it  is  assumed  that  there  is  no  y-dlrectlon  varia¬ 
tion  of  the  boundary  undulations). 

-  exp  i(cx  -  uQz)  (IV.1) 

where  oq  and  are  the  projections  of  the  incident  wave  number  upon  the  x-  and 
z-axes  respectively,  i.e.,  ao»kosin0o  and  ^o»kocos0o(ko*»co/co). 

According  to  Rayleigh,^  this  incident  wave  excites  a  denumerably  infinite 
set  of  plane  reflected  waves  <Pr  and  of  plane  transmitted  waves  (p^,  which  can  be 
written 

00 

<Pr  ■  /  An  exp{i[(co  +  ga)x  +  pnz])  (IV.2) 

with  n2  -  k2  -  (cq  +  gn)2 

*The  harmonic  time  factor  exp(-ln)t)  will,  as  usual,  be  suppressed. 
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(IV. 3) 
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Bn  exp{l[(«0  +  gn)x  -  X,,*)) 


with  V  -  +  gn)2 


(IV.  4) 


(IV.5) 


where  g  is  proportional  to  the  reciprocal  of  the  wavelength  of  the  boundary 
undulations  (see  Fig.  3)«  Kef.  16  contains  an  excellent  physical  explanation 
for  writing  the  reflected  and  transmitted  waves  in  this  fashion.  The  restriction 
Eq.  IV. 3  on  the  direction  cosines  of  the  reflected  rays,  nn,  results  from  the 
fact  that  Eq.  IV. 2  must  satisfy  the  wave  equation  in  the  incident  fluid*  the  re¬ 
striction  IV. 5  arises  in  an  analogous  fashion. 

The  ray  defined  by  n*0  is  the  specularly  reflected  one,  l.e.,  the  one  for 
which  the  angle  of  reflection  equals  the  angle  of  incidence  and  the  only  one  which 
would  be  present  if  the  boundary  were  plane.  The  rays  defined  by  n  >  0  are  some¬ 
times  called  forward  scattered,  while  those  for  n  <  0,  backward  scattered.  If 
2  2 

(o+gn)  is  greater  than  kQ,  for  some  n,  then  the  corresponding  nn  is  imaginary  and 
hence  that  ray  and  all  rays  defined  by  larger  values  of  n  are  exponentially  atten¬ 
uated  with  distance  from  the  boundary  and  hence  contribute  nothing  to  the  sound 
field  at  large  distances  from  the  boundary;  these  attenuated  waves  are  called 
inhomogeneous  or  evanescent  waves.  Note  that  if  ^  is  less  than  (l.e.,  g  >  k) 
all  rays  except  n«0  (and  possibly  n-1,  depending  on  0q)  are  inhomogeneous,  i.e., 
the  rough  boundary  reflects  Just  as  a  plane  boundary. 

The  unknown  coefficients  A  and  B  are  determined  from  the  two  boundary  con- 

n  n 

ditions  which  are  continuity  of  pressure  across  the  boundary  and  continuity  of 
normal  velocity. 

po<P  »  (IV. 6) 


dtp/dn  «  dcp^/dn 


(IV. 7) 


where  «p»9i+9r  and  n  indicates  the  normal  direction  to  the  boundary.  The  equation 
of  the  boundary  is  z«((x)  consequently  Eq.  IV. 7  can  be  written 


d(9  -  <PX) 
5z 


d(<P  -  <PX) 
5x 


0 


(IV. 8) 


Q 

At  this  point  in  the  analysis,  we  must  make  the  basic  approximation  Rayleigh 
did.  It  consists  in  saying  that  the  slope  of  the  boundary  is  very  small,  i.e., 
d£/dx  «  0. 
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For  £(x)=a  cos(2rtx/\b),  this  means  that  wO  or  the  amplitude  of  the 

undulations  must  be  very  much  smaller  than  their  wavelength.  With  this  approxi¬ 
mation  and  the  definitions  of  the  potentials  in  Eqs.  IV.  1,  IV. 2,  and  IV. 4,  the 
boundary  conditions  become  (the  limits  of  summation  which  are  always  -®  to  +® 
are  omitted). 


BnexP  i 1 1  (°0+gn)x-xMZ  1 ) 


=  "o 


+  y^Anexp{i[(aQ+gn)x+|i, 


n  "'Z=C(x)  pl 

4. 


T“  |exp{i(aox-uoz)} 


C(x) 


(IV. 9) 


■i^exp  (i (a^x-p^z ) }+  i  ^An^nexp{i[(a+gn)x+|.t„z]) 


o  o 


n 


+  i  .  B  Y  exp{i[(a  +gn)x-x  z]} 
n*n  o  °  '  n 


z-C(x) 


=  0 


(IV. 10) 


At  this  point  in  the  analysis  we  must  introduce  another  approximation  which 
naturally  restricts  the  applicability  of  the  results  but  which  makes  the  solution 
of  Eqs.  IV. 9  and  IV. 10  tractable.  This  approximation  consists  of  restricting  our 
attention  to  high  frequencies  such  that  \q  «  It  is  then  true  that 
gn  «  a  <  k.  We  also  assume  that  the  wavelength  in  the  transmitted  medium  is 
much  less  than  so  that  gn  «  k^.  It  is  also  evident  that  for  these  inequali¬ 
ties  to  hold  we  must  restrict  our  attention  to  rather  small  values  of  n,  i.e., 
to  the  first  one  or  two  rays  on  either  side  of  the  principal  (n=0)  ray.  However, 
as  noted  earlier,  the  rays  defined  by  large  n  are  quite  likely  to  be  inhomogeneous 
and  thus  contribute  nothing  to  the  sound  field  at  a  large  distance  from  the 
b  oundary. 

With  the  approximations  just  described,  we  can  write 


X 


PS 


n 


n 


\f 2  21 

vk  -a  =  k  cose 
0  0  0  0 

^k 2 -a2'  =  k  cose 
1  o  1  1 


(IV. 11a) 
(IV.  lib) 


Snell's  law,  (sine^)/c^=(sineo)/co  was  used  to  obtain  Eq.  IV. lib.  We  can  now 
bring  and  ^  out  from  under  the  summation  signs  in  Eq.  IV.  10. 

In  this  report,  we  are  primarily  concerned  with  determining  what  is  reflected 
from  the  boundary  and  hence  we  will  only  solve  for  the  coefficients  Ar,  although 
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there  is  no  difficulty  in  solving  for  the  transmitted  rays  as  well.*  We  use 
Eq.  IV. 9  to  eliminate  the  coefficients  Bq  from  Eq.  IV.  10  which,  after  a  little 
algebra,  can  be  written 


where 


Vzo 


Zl+Z0 


exp[-i  2kQ£(x)  cos0 J  +  ^  Aq  expfignx]  *  0 


plCl 

COS0, 


and 


P  c 
o  o 

COS0 


(IV.  12) 


(IV. 13) 


The  first  factor  in  the  first  term  of  Eq.  IV.  12  is  of  course  the  plane  boundary 
reflection  factor  which  we  denote  by  R.  Employing  the  standard  identity 


exp[-i  a  cos(gx)] 


(-i)n  Jn(a)  exp[ignx] 


(IV.  14) 


in  Eq.  IV. 12,  where  J  is  the  Bessel  function  of  the  first  kind  of  order  n,  we 

n 

finally  obtain  this  result. 

R  Y'  (-i)n  J  (2k  a  cos0  )  exp(ignx]  ■  ^  A  exp[ig.jc].  (IV.15) 

L _  no  o  c _ *  n 

It  is  now  a  trivial  task  to  obtain  the  coefficients  Ar  by  equating  terms  with 
the  same  x-dependence  [note  that  J  n=(-l)  JnJ. 

A  =  RJ  (2k  a  cos0  )  A0  =  -RJ  (2k  a  cos0  ) 

o  ox  00  2  2  o  o 

A,  =  -iRJ^ (2k  a  cos0  )  A  0  *  -Rj(2k  a  cos0  ) 

1  1'  O  O  -2  2  O  O 

A  ,  =  -iRd, (2k  a  cos0  )  etc. 

-1  1  o  0 

These  results  are  very  similar  to  those  obtained  by  Rayleigh  for  the  case 

9 

of  normal  incidence;  in  fact  Rayleigh,  before  leaving  the  problem  of  scattering 
from  a  rough  boundary,  gives  the  above  result  for  the  amplitude  of  the  specularly 
reflected  ray  (i.e.,  n=0).  It  is  noted  that  the  amplitude  of  the  corresponding 
rays  on  either  side  of  the  specular  direction  are  precisely  equal,  i.e., 

A^=A  y  The  results  Rayleigh  quotes  for  normal  incidence  are  presumably  the  total 


*In  Sec.  B  for  illustrative  purposes  we  actually  solve  for  other  than  the 
reflected  waves. 
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amplitude  of  the  pair  of  corresponding  rays  since  his  is  twice  the  result  given 
here  and  he  quotes  no  result  for  A  For  large  argument,  all  Bessel  functions 
look  like  damped  sinusoids;  therefore,  with  increasing  frequency,  less  and  less 
sound  is  reflected  in  the  specular  and  near-specular  directions.  This  is  one  of 
the  basic  features  of  the  transient  response.  Furthermore  the  most  damaging  com¬ 
ponent  of  the  pressure  wave  is  associated  with  the  specular  reflection. 

We  turn  now  to  the  consideration  of  scattering  of  a  plane  wave  incident  from 
a  fluid  onto  other  media.  It  is  stressed  that  throughout  the  remainder  of  the  re¬ 
port,  the  two  basic  approximations  already  encountered,  i.e.,  that  the  slope  of 
the  boundary  is  small  and  that  the  wavelengths  of  all  reflected  or  refracted  com- 
pressional  or  shear  waves  are  small  compared  to  the  wavelength  of  the  boundary, 
will  be  employed.  Because  the  analyses  are  very  similar  to  the  preceding,  their 
presentation  will  be  quite  brief. 


B.  Scattering  from  a  Sinusoidal  Boundary  between  a  Semi- Infinite  Fluid  Medium 
and  a  Semi-Infinite  Elastic  Medium 

The  geometry  of  the  problem  is  essentially  that  of  Fig.  3  except  that  we  now 
assume  the  upper  medium  to  have  a  shear  wave  velocity  b^  and,  correspondingly,  a 
set  of  shear  waves  defined  by  the  potential  ^ 

00 

*  Cn  exp{i[(aQ+gn)x-Jnz]}  (IV.16) 

n=  -  ® 


2  2  (  A  x2 

»„  ”  H1  -  (Vgn) 


(IV. 17) 


where 

H1  =  “^l  * 

The  boundary  conditions  are  continuity  of  normal  displacement,  continuity 
of  normal  stress  and  vanishing  of  the  tangential  stress.  In  terms  of  the  potentials 
cp  (=<Pi+9r),  «PX  and  ijr^,  these  become^  (assuming  d£/dx  «  0) 


dco  ^1  b'l 

•l  +  sr  +  3r  =  0 

2  2 

p  0  -d  d  (If, 


at  z=£(x) 


(IV. 18) 


82cpx 

8xoy 


a2* 

az2 


1  d  *1 
"  +  =  0 
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The  assumption  that  all  wavelengths  are  small  compared  to  that  of  the  boundary 
allows  us  to  write,  as  before. 


p  ms  k  cos0 
n  o  o 


cos0^ 

v  «  k,  cosy, 
n  1  T1 


(!V.19a,b,c) 


The  angles  0^  and  y^  are  defined  by  Snell's  law 

sin0  sin0,  siny.. 
oil 


(IV.20) 


where  y  is  the  angle  that  the  principal  (n=0)  ray  of  the  set  of  transmitted  shear 
waves  makes  with  the  z-direction. 

For  convenience,  let  us  introduce  the  notation 


A  =  ^An  exp[i(gnx  +  nnC)3 
B  =  ^  Br  exp[i(gnx  -  y^Q)] 
C«£cn  exp[i(gnx  -  unC)] 


(lV.21a,b,c) 


Whence  the  boundary  conditions  Eq.  IV.  18  become,  after  some  algebraic  manipulation 


-  COS0  - -  COS0. 

o  c^  1 


.  c  sin  2y. 

—  -  A  -cos0  exp(-ik  Qcosd  ) 

2  b1  cosy1  o  o*  o 


-  cos  2y. 

pi  1 


-  sin  2y^  I  B  ■  — -  exp(-ikoCcos0Q) 


(IV.22) 


b1  sin  2y1 
cq  cosy^ 


COS0 


^  cos  2y1 


Solving  for  A,  we  find 


where 


A  ■  exp(-ik  £cos0  ) 
o  o 


pibi 


t  cosy. 


"2  2 

Z^cos  2y^  +  Z^sin  2y^  -  Zq 

p  2 

Z^cos  2y^  +  Ztsin  2y^  +  Zq 


(iv. 23) 


(IV.211 


while  Z  and  Z,  were  defined  in  Eq.  IV.  13. 
o  1 
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The  term  in  brackets  is,  again,  simply  the  plane  boundary  reflection  factor 
(denoted  V  in  Brekhovskikh^).  Therefore  we  write 

Anexp(ignx)  =  V  exp[-i(2kQa  cos0o)cos(gx) ]  (IV.25) 

which,  with  the  identity  IV. 4,  becomes 


QD 

I 

n=  -® 


Anexp(ignx) 


=  V 


£  <->\ 


(2kQa  cos0q)  exp(ignx) 


n=  -« 


whose  solutions  are 


A  =  V  J  (2k  a  cos 0  ) 
0  0  o 


A, ,  =  -i  V  J, (2k  a  cos0  ) 
±1  1  o  0 

A. „  =  -V  J_(2k  a  cos0  ) 
x 2 !  2  o  o 


(IV. 25) 


(IV.27) 


Note  that  if  we  solve  Eq.  IV.22  for  B,  we  obtain 

y  Bnexp(ignx)  =  W  exp{-i[(kocos0o  -  k-^cos©^  a  cosgx]}  (IV. 28) 

where  W  is  the  amplitude  of  the  transmitted  compress ional  wave  in  the  case  of 
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a  plane  boundary.  Eq.  IV. 28  has  the  solutions 


Bn  =  V  (-i)nJn[(k0  cos0Q  -  k^cos0^  )]  •  (IV.29) 

If  we  solve  Eq.  IV.22  for  C,  we  obtain 

i 

1  CRexp(ignx)  =  P  exp{-i[(kQcos0o  -  H^cosy^Ja  cosgx])  (IV.30) 

where  P  is  the  amplitude  of  the  transmitted  shear  wave  in  the  case  of  a  plane 
19 

boundary.  Eq.  IV. 30  has  the  solutions 


CR  =  P  (-i)nJn[(kocos0Q  -  H1cosy1)a  ]  (IV.31) 

C.  Scattering  from  the  Sinusoidal  Surface  of  a  Thick  Fluid  Layer  Sandwiched 
between  Two  Different  Semi -Infinite  Fluid  Media 

The  geometry  of  the  problem  is  illustrated  in  Fig.  4.  In  addition  to  the  in¬ 
cident  wave  9.  ,  the  reflected  wave  in  the  incident  medium  9  and  the  transmitted 
j.  r 

wave  in  the  second  medium  9j_>  must  allow  for  a  reflected  wave  from  the  plane 
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/ 


surface  of  the  intermediate  layer  (the  surface  at  z=  -h)  which  we  denote  cp2  and 
for  a  wave  transmitted  into  the  third  fluid  media,  which  we  denote  cp^ .  We  write 
these  two  additional  potentials  in  this  fashion 


qp 

*  Dnesp(l[co+gn)x  +  x^ll 


(IV.32) 


2  2  2 
where  ^  =  ^-(a^gn) 


(IV.33) 


00 

cpA  =  \  B  1 

A  B 

n=-® 


exp{i[(aQ+gn)x  -  t^z]} 


(IV. 34) 


where  \  m  *2  ~  (°0+8n)2 


(IV.35) 


It  is  evident  that  even  though  we  have  a  plane  boundary,  these  two  wares  must 
be  written  in  the  same  fashion  as  the  infinite  set  of  waves  associated  with  a 
sinusoidal  boundary  to  be  able  to  match  the  boundary  conditions;  for  example,  if 
the  plane  boundary  of  the  intermediate  fluid  were  a  free  boundary  the  potentials 
cp^  and  <P2  would  obviously  have  to  cancel  each  other  and,  given  tp^,  (Eq.  IV. 4), 
this  could  only  be  possible  if  has  the  form  defined  above. 

The  boundary  conditions  are  again  continuity  of  pressure  and  normal  velocity 


become 


£(x)  and  z* 

-h.  Under  the 

P0(«Pl  +  9r) 

«  P1('P1  +  <P2) 

+  <Pr) 

+  «p2) 

dz 

dz 

Pl({pl  +  *2* 

"  Pa’t 

dfa-L  +  <P2) 

dz2 

'  5T 

at  z »C(x) 


(lV.36a,b) 


at  z=  -h 


(lV.37a,b) 


As  usual,  the  assumption  that  all  wavelengths  are  small  compared  to  that  of  the 
boundary  allows  us  to  approximate  the  direction  cosines  as  in  Eq.  IV.  19a, b  and 
in  addition  to  write 


\  «  k2cos02 


(IV. 38) 
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where  d 2  is  defined  by  Snell's  law 


sind  sind.  sind_ 

o  1  2 


(IV. 3?) 


Employing  the  shorthand  notation  introduced  in  Eq.  IV. 21  (where  D  and  E  are  the 
right-hand  sides  of  Eqs.  IV. 32  and  IV.33>  less  the  factor  exp(ipz),  the  four 
boundary  conditions  can  be  cast  in  the  form 

Ple-lklCcosaI  PiettlCcosSl  0 

k  cose  e^oCcoseo  -IkiCcosex  ikiCcosej. 

o  o  11  11 


^ik^hcosd^ 


4-ik^hcosdjL 


•P2e 


ik2hcosd2 


k1cosele-Utlhcos9l 


-k1cose1e-Utlhcos9l  -k2cose2elk2hco89d  E 


ik  Ccosd 
p  e  o  o 
Ko 

,  .  -ik  Ccosd 

k  cosd  e  o  o 
o  o 


(iv.4o) 


Solving  for  A,  the  amplitudes  of  the  set  of  reflected  waves  in  the  incident  medium 


we  obtain 


JO 

2> 


eignx=e-iZkoCcos0 


n=  -® 


where 


z  z  z 

(-1+  2^)cos[k1(htC)cosd1]-i(^  -  ^)sin[k1(htC)cosd1J 

2  Z  Z 

(1+  ^)cos[k1(h+C)cosd1]-i(^:  +  ^)sin[k1(htC)cosd1] 


(IV. 41) 


Z  -  ^ 
2  -  cosd,. 


(IV. 42) 


while  Z^  and  Z^  were  defined  in  Eq.  IV.13* 

Now  if  Q  is  negligible  compared  to  h  (i.e.,  if  the  amplitude  of  the  boundary 
undulations  are  small  compared  to  the  thickness  of  the  intermediate  layer)  so  that 
we  may  discard  it  in  the  arguments  of  the  sines  and  cosines  in  the  factor  in  brackets 
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in  Eq.  IV.  41,  then  that  factor  reduces  to  the  standard  plane -boundary  reflection 
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factor  for  the  three -fluid  layer  problem.  We  denote  the  factor  by  0. 

Employing  the  Bessel  function  identity  Eq.  IV.  14,  Eq.  IV.  41  becomes 


00 


n=-® 


G 


£ 

Qs  .00 


(-i)neignx  J 


(2k  a 
nx  o 


COS0q) 


(IV.  43) 


which  has  solutions 


A  =  G  (-i)n  J  (2k  a  cos0  ) 
n  n  o  o 


(IV.  44) 


If  one  were  to  solve  for  the  amplitudes  of  the  various  other 
B,  D,  or  E,  the  results  would  quite  evidently  be  the  appropriate 
transmission  factor  from  the  plane  boundary  case  multiplied  by  a 
whose  argument  involves  the  difference  of  wave  numbers  as  in  Eq. 


sets  of  waves, 
reflection  or 
Bessel  function 
IV. 29  and  IV. 31. 


D.  Scattering  from  the  Sinusoidal  Surface  of  a  Thick  Elastic  Layer  Sandwiched 
between  Two  Different  Semi-Infinite  Fluid  Media 

The  geometry  of  the  problem  is  the  same  as  in  Fig.  4  except  we  now  assume  the 
intermediate  layer  to  have  the  shear  wave  velocity  b1  and  correspondingly,  two 
sets  of  shear  waves  -  one  set  that  originates  at  the  sinusoidal  boundary  and  the 
other  that  represents  reflections  from  the  plane  boundary. 

It  is  evident  from  the  preceding  work  that  under  the  three  basic  assumptions 
made  in  this  analysis,  viz.,  negligible  boundary  slope,  all  wavelengths  much  less 
than  the  boundary  wavelength  and  layer  thickness  much  greater  than  amplitude  of 
boundary  undulations,  the  reflection  or  transmission  factors  for  a  sinusoidal 
boundary  are  simply  the  corresponding  factor  for  a  plane  boundary  multiplied  by 
a  Bessel  function  whose  argument  involves  the  wave  numbers  in  the  appropriate 
media,  the  amplitude  of  the  boundary,  and  the  angle  of  incidence.  We  shall  there¬ 
fore  spare  the  reader's  time  and  simply  quote  the  answer  for  the  reflection  co¬ 
efficients  Ar  of  the  reflected  potential  (Eq.  IV.2), 

A  =  T(-i)n  J  (2k  a  cos0  ) 
n  x  o  0  o 

where  T  is  given  in  Ref.  21. 
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E*  Spattering  from  a  Sinusoidal  Boundary  between  Two  Semi-Infinite  Fluid 
Media  -  Slope  of  Boundary  Not  Negligible 


Let  us  return  to  the  problem  considered  in  Sec.  A  and  relax  the  condition 
on  the  slope  of  the  boundary,  i.e.,  let  us  try  to  solve  the  problem  using  the 
exact  form  of  the  second  boundary  condition  (Eq.  IV. 8).  For  £=a  cosgx,  this 
second  boundary  condition  can  be  written 


-H  0-^7 A  ,i  ei(gnx+Mn')  +  Y  B  x  = 

o  L — >  n'  n  l~  nAn 

^afi[elg>c.e-1e>c][7Bn(-o+en)eifeIa-xnh  -ae_i^o^-  7An(a0+8n)ei^8“t+llnC  )  ^ 


Employing  the  approximation  that  all  wavelengths  are  much  smaller  than  that  of 

the  boundary  so  that  u  cos 0  ,  \  »k.,  cos0,  and  a  +gn%k  sin0  (note  that  because 

no  o  n  l  1  o  o  o 

of  this  last  approximation,  the  results  below  will  not  be  strictly  valid  for  nor¬ 
mal  incidence),  eliminating  the  B  from  Eq.  IV. 45  by  the  use  of  Eq.  IV. 9,  and 
finally  employing  the  Bessel  function  identity  Eq.  IV. 14,  we  obtain  this  relation¬ 
ship  for  the  amplitudes  of  the  reflected  rays 


<wD-i,n  V2koa  case>a)  ei3nx + (wL  vis“  ■ 


p[Zlta^o(^  .  [(-!>"  Jn(2Koo  cos9o)+An]eignX 


(IV.U6) 


-ia 

o 


If  we  now  equate  terms  with  the  same  x-dependence  on  both  sides  of  the  equation, 
we  generate  an  infinite  set  of  coupled  equations  each  of  which  involves  three  of  the 
unknown  A^.  For  example,  gathering  the  terms  that  have  no  x-dependence  and  that 
have  x-dependence  exp(±igx)  results  in  these  three  equations 


A0  ■  M<“>  +  1-A1+A-!J 

A1  •  +  [J0(a)w2(a)+A0-A2] 


(lV.47a,b,c) 


where 


A-l  "  -1RJl(a)  +  z^ 


1+Z0 


a  =  2k~a  cos0  and  3  = 
0  o 


[-J(a)-J2(a)-A0+A_2] 
(-iag/2)^  tan0o(po/p1-l)]. 
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t 


In  general,  the  nth  coefficient  A  is  coupled  to  both  A  .  and  A .  _ .  We  solve 

n  n-±  Dri. 

these  by  iteration,  i.e.,  we  use  the  results  obtained  in  Sec.  A  (now  denoted  A 
in  the  right-hand  side  of  these  equations,  to  generate  a  first  order  correction 
(denoted  A^  which  includes  the  slope  of  the  boundary.  Noting  that  A^*A^, 


(0) 


we  obtain 


aW  =  Rj(a) 

fi)  (22. ) 

*£'  =  -lW^a)  +  0 - [JQ(a)  +  J2(a)] 


-  -iBjjfa)  -  S 


(zl+z o)‘ 

(2\) 


[J.(a)  +  J0(a)] 


2  1  o 


j-imploying  the  identity 

Jn-l^a)  +  Jn+l^a)  =  2nJn^a)/a  ' 
the  last  two  equations  can  be  written 

a±i>  ■  r£se 


[^(Pq/Pj^)] 


bCOBOo  Ci-(z0/zir] 


_  tand 
2,  o 


(lV.48a,b,c) 


(IV.  1*9) 


(IV.50) 


In  general 


(-i)“EJ  (2k  a  cose  ){1V  °fi  (pc/pl^  tan0 

n  o  0  koco8®o  [1-(Z 0 


(IV. 51) 


The  correction  term  is  seen  to  be  independent  of  a,  the  amplitude  of  the  boundary 
undulations.  Also  the  correction  term  has  destroyed  the  symmetry  of  the 
amplitude  coefficients  so  that  A_n  no  longer  equals  An> 

V.  THE  REFLECTION  OF  PLANE  WAVES  FROM  A  LIQUID -SATURATED  POROUS  ELASTIC  HALF-SPACE 

In  an  earlier  report'  formulas  were  derived  for  the  reflection  coefficient 

for  a  plane  harmonic  wave  incident  on  a  porous  fluid  saturated  elastic  solid  using 
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field  equations  derived  by  Biot.  The  reflection  coefficient  so  formulated  de¬ 
pends  on  seven  parameters  characterizing  the  liquid-saturated  porous  elastic  solid 
as  well  as  the  density  and  sound  speed  of  the  fluid  in  which  the  plane  wave 
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propagates.  The  seven  parameters  characterizing  the  porous  elastic  solid  are: 

1.  p i  ,  the  density  of  the  pure  compacted  solid. 

2.  p^.,  the  density  of  the  fluid  which  saturates  the  solid. 

3.  Cu,  the  unjacketed  compressibility  of  the  liquid  saturated  elastic  solid. 

4.  Cj,  the  jacketed  congress ibility. 

5»  N,  the  shear  modulus  of  the  material. 

6.  "the  structure  factor. 

7»  P,  the  porosity. 

The  physical  significance  of  the  above  parameters  are  described  in  the  report 
by  Eichler  and  Rattayya.1 

Certain  of  the  above  parameters  can  be  determined  from  experiments  while  the 

others  must  be  estimated.  For  example ,  the  only  measurements  made  of  un Jacketed 
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and  Jacketed  compress ibility  were  those  made  by  Fatt.  However  these  were  made 
for  a  kerosene  saturated  sandstone  under  static  conditions  in  the  laboratory  and 
in  any  case  would  not  be  applicable  to  a  theory  which  attempts  to  model  the  ocean 
bottom.  A  similar  situation  exists  with  regard  to  x>  the  structure  factor;  a 

2k 

factor  x  approximately  equal  to  4.3  has  been  measured  for  an  air- filled  sand. 

Very  few  measurements  of  the  shear  wave  velocity  in  ocean  bottom  sediments  are 
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available.  In  fact  the  only  oner,  we  are  aware  of  are  those  reported  by  Laughton. 

With  the  above  limitations  in  mind,  a  numerical  procedure  to  determine  the 
reflection  coefficient  as  a  function  of  the  grazing  angle  has  been  developed  and 
programmed  for  a  digital  computer.  Assumed  values  for  the  structure  factor  x 
and  the  unjacketed  compressibility  together  with  experimentally  obtained  values 
of  the  porosity  0,  saturating  fluid  density  pf  (in  all  cases  the  saturating  fluid 
density  is  that  of  sea  water)  the  bulk  density  p  (mass  per  unit  volume  of  the  com¬ 
bined  solid  and  liquid),  the  dilatational  wave  speed  V  ,  and  the  shear  wave  speed 

V  ,  are  entered  into  the  program.  The  program  then  proceeds  to  calculate  the 
s 

quantities  necessary  for  the  determination  of  the  reflection  coefficient. 

In  each  case  the  calculations  were  performed  for  a  relatively  large  structure 
factor  x“H  sod  "two  smaller  structure  factors  x=l,  and  x=2*  The  data  used  was 
that  of  Laughton  on  Globigerina  Ooze. 

In  both  Figs.  5  and  6  we  show  plots  of  the  bottom  loss  per  reflection  versus 
grazing  angle  for  a  given  set  of  parameters.  The  bottom  loss  per  reflection  is 
defined  by  [-20  log|R|]  where  |R  |  is  the  modulus  of  the  reflection  coefficient. 
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In  each  figure  we  show  two  plots,  one  for  a  value  of  Cu=.l  Cf  and  the  other  for 
Cu«.2  Cf,  where  is  the  compressibility  of  the  saturating  fluid.  Fig.  7  shows  the 
bottom  loss  per  reflection  for  a  given  set  of  parameters  for  two  different  values 
of  the  structure  factor  x=l>  and  X=2.  It  is  characteristic  of  each  of  these  plots 
that  there  is  some  angle  at  which  the  bottom  loss  is  a  maximum.  There  is  also 
some  angle  above  which  the  bottom  loss  is  relatively  insensitive  to  change  in  the 
grazing  angle.  Fig.  8  shows  a  plot  of  the  phase  angle  cp  of  the  reflection  coef¬ 
ficient  superin$osed  on  a  plot  of  bottom  loss.  From  this  figure  we  see  that  the 
angle  above  which  the  bottom  loss  is  fairly  constant  is  precisely  the  critical 
angle  9 c,  where  is  given  by  0c*cos~  (c0/^p)>  c0  the  sound  speed  in  the  fluid 
and  V  is  the  dilatational  velocity  in  the  solid  (i.e.,  the  angle  for  vhlch  the 
phase  angle  of  the  rjflection  coefficient  goes  to  zero). 

Since  we  have  neglected  viscosity  in  the  above  calculations  the  reflection 
coefficients  obtained  are  independent  of  frequency. 


APPENDIX  A:  Laplace  Transform  of  Time  Dependent  Green’s  Function  for  a  Wedge 
with  Pressure-Release  Walls 


Consider  the  problem  of  a  source  with  Dirac  delta  time  dependence  located  at 
the  point  S(ro,0Q)  within  a  wedge-shaped  region  0  <  G  <  a,  the  boundaries  of  which 
are  pressure-release  surfaces.  The  pressure  must  therefore  satisfy  the  equation: 


dj>  1  dg  ,  1_  _  1_  dfg  «(r-*,,)»(Me)t(t) 

dr2  r  "  '  r2  da2  c2  at2  ’  r 

and  the  boundary  conditions 
p(r,0)  -  0 
p(r,a)  «  0  . 


(A.l) 


(A.2) 


Taking  a  Laplace  Transform  in  time  yields  the  following  system  of  equations: 

ii.  .15  *1.  ii  l-  ‘fr-roW9-V 

Sr2  75?  r2*2'caP“ 

and 

p(r,0)  *=  0 
p(r,a)  =  0 

where  p(r,0)  is  the  Laplace  Transform  of  p(r,0,t). 


(A.3) 


(A.4) 
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The  angular  eigenfunctions  corresponding  to  the  boundary  conditions  (B. 4)  are 
V0)  =  An  sin  ^  (A-5) 

where 


Dn  =  ,  n=l,2,3  ••• 


(A.6) 


Assuming  that  the  transformed  pressure  p  can  be  expanded  in  the  form 


OP 

p(r,0;s)  -7  An(r;s)  sin 
n=l 


(A. 7) 


and  using  the  orthogonality  of  the  eigenfunctions  over  the  range  0  <  0  <  a  it 

is  found  that  the  functions  A  must  satisfy 

n 

d2A  ,  dA  2  v2  0  6(r-r  ) 

_ n  ,  1  n  r s  ,  n\  A  _  2  a  _ _ o' 

,2  r  dr  "  ^  2  2^  n  ”  ~  Cl  s  un0o  r 

dr  c  r 


(A.8) 


Equation  A.o  is  solved  using  Hankel  Transforms  and  the  result  is 


A  (r ,s)  =  -  sin  u  0  I  (-  r  <)  K  (fr  >) 
nx  ’  a  n  o  u  c  c 

n  n 


(A.9) 


where  r  <  and  r  >  are,  respectively,  the  smaller  or  larger  of  the  quantities 

r  and  r  . 
o 

Using  the  identidy 

sin  v  0  sin  v  6  =  ^  {cosv  (0-0  )  -  cosy  (0+0  )} 
no  n  2  n'  o'  nN  o 

the  solution  to  the  set  of  equations  A. 3  and  A. 4  can  be  written  as 


(A. 10) 


p(r,0;s) 


w 


e„ [c°sv„ (e-e„)-cosvr,(e+0„) ]  I,  (f  r  <)K  (£  r  >)  (A.ll) 

n  n 


2a/  *nl  n'~  ’o'  - n'~  o' '  "u_'c  '  u_'c 

fe) 


where  e  =1  and  e  =2  for  n  >  1. 
on  — 
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APPENDIX  B:  Asymptotic  Evaluation  of  I 


Consider  the  integral  I  where  I  is  of  the  form 

I  -  £  J f  Aft)  )  dS 

>QO 


(B.l) 


where  the  integral  is  taken  along  the  real  axis  in  the  £ -plane  and  the  branch  of 
is  that  for  which  Im  /k*  -?2  >  0  (see  Fig.  9a)* 

We  define  a  conformal  transformation  5=kQsin<p  from  the  § -plane  to  the  region 
-  ^  <  Recp  <  ~  of  the  <p-plane.  Then 


1  =  ±  J*  A(kosincp)  ei*0(xsin<f*zcoscp)  kQcoscp  dtp  ) 
c  ° 

=  J  A(k  sin«p)  eikoCOS^-0^  costp  dtp  j 

dX  °  O  ' 

c 


(B.2) 


where  x=r  sin0,  and  z=r  cos0,  and  where  the  contour  C  is  as  shown  in  Fig.  9b. 

This  integral  may  be  evaluated  asymptotically  for  the  case  where  kQr  »  1  by  the 
method  of  steepest  descents.  Letting 

F(<p)=  coscp(kQsin<p)  (B.3) 

we  have 

1  =  2k  J  F(*P)  e^o^06^-0)  dtp  *  (B.4) 

c 


The  location  of  the  saddle  point  is  found  from  the  condition 

[S^ikTCOSKHO]  =  0 


(B.5) 


so  that  the  saddle  point  (SP)  is  at  tp  =  0. 
through  the  saddle  point  is  given  by 

Im  {ikQr  cos(tp-e)  -  ikQr}  =  0 

i.e.,  Re[cos(tp-@)-l3  =  0  > 

or  cos(u-0)  cosh  v  =  1 

/ 

if  tp  =  u  +  iv. 


The  path  of  steepest  descent  (PSD) 


(B.6) 
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This  path  is  sketched  in  Fig. 9b.  We  deform  the  contour  C  into  the  PSD  and, 
since  the  integrand  decays  rapidly  away  from  the  saddle  point  on  this  contour,  we 
approximate  the  integral  by  its  value  along  a  small  section  of  the  path  around  the 
saddle  point.  The  leading  term  of  the  asymptotic  expansion  for  I  is  then  given  as 


I  ~ 


cose  A(kQsine)  e^k0r"*A) 


(B.7) 


provided  9  fa/ 2,  and  Arsine)  is  bounded  and  non-zero.  To  this  must  be  added 
the  residues  of  the  integrand  at  all  poles  which  lie  between  the  original  path  C 
and  the  steepest  descent  path  (a  possible  set  of  these  poles  are  shown  circled  in 
the  figure).  The  poles  are  given  by  the  zero3  of  the  transcendental  equation 


i  tan  y^h 


Vi 

Vo 


where  y^ 


(B.8) 


Since  A^(5)  is  an  even  function  of  y^,  there  are  no  branch  cuts  necessary 
to  make  y^  single -valued,  so  long  as  h  is  finite.  As  h  -»  ®  the  poles  become  more 
and  more  numerous  along  a  line  which  in  the  limit  becomes  a  branch  line. 
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Fig.  2.  Geometry  of  the  problem  associated  with  a  plane  wave 
<P^  incident  on  an  arbitrary  rough  interface  z=«f(x). 
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Fig*  3*  Geometry  of  the  problem  associated  with  a  plane  sound 
wave  obliquely  incident  onto  the  sinusoidally  rough 
boundary  between  two  semi-infinite  fluid  media. 
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fig.  4.  Geometry  of  the  problem  associated  with  a  plane  sound 
wave  obliquely  incident  onto  the  sinusoidal  boundary 
of  a  thick  fluid  layer  sandwiched  between  two  different 
semi-infinite  fluids. 
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£  Plane: 


4>  Plane: 


Fig.  9.  (a)  Integration  path  for  evaluation  of  integral  in  Eq.  B.l  in  the 

complex  5 -plane.  (b)  Integration  path  in  the  complex  cp-plane  for 
evaluation  of  integral  in  Eq.  B.2,  also  showing  steepest  descent 
path  as  well  as  possible  pole  locations. 
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